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ABSTRACT 

We study the cosmological evolution of BL Lac objects by applying the analysis to a new 

sample of 55 objects presented for the first time in this paper. This sample has been selected from 
the 239 sources with the brightest X-ray flux (>4xl0~^^ erg s^"'^ cm~^) and relatively bright optical 
counterpart (B<20.5) among the ^^1600 objects included in the REX survey. The uniform distribution 
of the Ve/Va values found in the sample suggests that BL Lac objects are not affected by a strong 
cosmological evolution in contrast with the behavior observed in the emission line AGNs. The VelVa 
analysis applied to the subsample of the High energy peaked BL Lacs (HBL) does not reveal a significant 
departure from a uniform distribution. This result suggests either that the cosmological evolution is less 
extreme than that previously found in other samples or that it is carried out only by a minority of 
objects, namely by the most extreme HBLs. 

Subject headings: surveys - galaxies: active - quasar: general - BL Lacertae objects: general 



, 1 Introduction 

■ Among the Active Galactic Nuclei (AGN), BL Lac ob- 
jects represent a minority characterized by an almost fea- 
iiureless optical spectrum and by an intense "blazar" activ- 
^ity (variability, flat radio spectrum, polarization). Due to 
their scarcity, the study of BL Lacs is not simple and many 
important issues, like the cosmological evolution, are not 
yet well understood. The first attempts to estimate the 
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.evolutionary behavior of BL Lac objects were based on 
'the Einstein Medium Sensitivity Survey (MSS, Maccac- 
|aro et aL 1982) and its extension (EMSS, Gioia et al. 
'1990). Evidence for negative evolution (objects less nu- 
merous/luminous in the past than now) was found (Mac- 
cacaro et al. 1984; Wolter et al. 1991; Morris et al. 1991; 
^Wolter et al. 1994; Rector et al. 2000). At the same 
'time, the analysis of the 1 Jy radio selected sample of 
'BL Lacs found a mildly positive evolution (Stickel et al. 
1991; Rector & Stocke 2001) at odds with the X-ray re- 
sults. However, all these conclusions were based on small 
samples containing 20-40 objects at most. 

The recent availability of large radio and X-ray sur- 
veys boosted the attempts at selecting large samples of 
BL Lacs. Different techniques are used in order to increase 
the efficiency of the selection process: a very efficient tech- 
nique turned out to be the combination of radio and X-ray 
information thanks to the fact that BL Lacs are both X- 
ray and radio loud sources as firstly pointed out by Stocke 
et al. (1991). This technique reaches efRciencies of I5%- 
30% (e.g. Wolter et al. 1997) depending on the flux limits 
involved, while in a purely X-ray or radio survey the per- 
centage of BL Lacs hardly reaches 5% (e.g. Stocke et al. 
1991). A radio survey of flat spectrum sources, like the 
I Jy survey, has an efficiency of about 12% (Stickel et al. 
1991) although deeper flat spectrum radio surveys include 
higher percentages of BL Lacs (up to 20% in the 200 mJy 



sample, Marcha et al. 1996 and in the CLASS Blazar Sur- 
vey, Marcha et al. 2001, Caccianiga et al. 2001). 

Thus, surveys based on the cross-correlation of radio and 
X-ray catalogs have been initiated (the DXRBS survey, 
Perlman et al. 1998, the RGB survey, Laurent-Muehleisen 
et al. 1998, the REX survey, Maccacaro et al. 1998; Cac- 
cianiga et al. 1999, hereafter Paper I). So far, the gen- 
eral picture is still quite unclear. Some surveys (Bade et 
al. 1998; Giommi, Menna & Padovani 1999; Giommi et 
al. 2001) find negative evolution but only for a restricted 
sub-set of the BL Lac objects, namely the ones with a syn- 
chrotron emission peaked at high energies (the High energy 
peaked BL Lacs, HBL). A continuous trend, ranging from 
slightly positive evolution for the Low energy peaked BL 
Lacs (LBL) to strong negative evolution for HBL has been 
found also by Rector & Stocke (2001) after the re-analysis 
of the 1 Jy and EMSS samples. At the same time, prelim- 
inary results from other surveys, like the DXRBS and the 
REX survey, are not finding any cosmological evolution, 
neither negative nor positive (Padovani 2001; Caccianiga 
et al. 2001). 

In this paper we present a study of the cosmological 
evolution of BL Lacs based on an almost (91%) complete 
sample of 55 BL Lacs created by imposing a relatively 
bright X-ray and optical cut on the REX survey. So far, 
this is the largest sample with a spectroscopic classifica- 
tion that is used to estimate the cosmological evolution 
of BL Lacs. The sample is presented in Section 2 while 
the statistical analysis is described in Section 3. In Sec- 
tion 4 we compare the result with those reported in the 
literature. Finally, the discussion and the conclusion are 
presented in Section 5. 
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2 The X-ray bright REXs sample 

The REX survey is the resuh of a positional cross- 
correlation between the NRAO VLA Sky Survey (NVSS, 
Condon et al. 1998) at 1.4 GHz and an X-ray catalogue 
of about 17,000 serendipitous sources detected in 1202 
pointed ROSAT PSPC fields (see Paper I for details). The 
good accuracy of the NVSS (VLA) positions is used to pin- 
point the optical counterpart. The spectroscopic identifi- 
cation of the sample is in progress. So far, about 40% 
of the sample has been spectroscopically identified and 
some of the newly discovered objects have been already 
presented (Wolter, Ruscica & Caccianiga 1998; Caccian- 
iga et al. 2000). The identification strategy that we have 
always adopted is to observe X-ray bright sources first, 
so as to obtain a complete subsample (of increasing size) 
with a very high fraction of identified sources usable for 
statistical studies. At present, at the X-ray flux limit of 
4xl0~^^ erg s~^ cm~^ and considering only the objects 
with a blue magnitude brighter than 20.5"*°^ the sam- 
ple is almost completely identified (identification rate of 
91%). We have thus analyzed this complete sample, that 
we call X-ray Bright REX (XB-REX), with the main goal 
of studying the statistical properties of the BL Lac objects. 
The incidence of the relatively bright X-ray cut-off on the 
final results is discussed in Section 4. Despite the fact 
that the REX survey has been created by using pointed 
PSPC images characterized by different exposure times 
and, thus, different flux limits, at the bright X-ray flux 
limit considered for the XB-REX sample the sky-coverage 
function is flat, i.e. the almost totality (99%) of the sur- 
veyed area reaches the X-ray flux limit of 4x 10~^^ erg s~^ 
cm^^. This greatly simplifies the analysis of the sample 
which is characterized by a single X-ray, radio and optical 
flux limit. 

In total, the XB-REX sample contains 239 objects. The 
magnitudes are preferentially taken from the APM (Auto- 
matic Plate Measuring^) database (O magnitude) or from 
the LEDA catalogue (Paturel et al. 1996). For about 10% 
of objects, for which no blue magnitude was available, we 
have estimated it from the APM E magnitude^ or used 
the B magnitude found in the literatTire. In case of the BL 
Lac objects presented in this paper, the great majority of 
the magnitudes are derived from APM except for 4 objects 
where the magnitudes are taken from the literature. 

In the XB-REX sample 91% of the objects (218 out of 
239) are already identified, either from the literature or 
from our own spectroscopy. In total, 55 BL Lac objects 
(25% of the identified objects), 95 emission line AGNs (EL 
AGN) and 66 objects classified as "galaxies" are found. 
The complete sample of 55 BL Lacs is reported in Ta- 
ble 2. In addition, two objects that are classified as stars 
are probably spurious optical counterparts and should be 
rejected from the sample (the real counterpart is likely 
to be fainter than the optical limit used for the XB-REX 
sample). 



The classification of the different classes follows the rules 
described in Paper I. In case of BL Lac objects, we have 
used the "expanded" definition based on the Ca II contrast 
at 4000A (A) which must be below 40% instead of 25% as 
firstly proposed by Stocke et al. (1989). This modification 
is justified by the evidence that the limit of 25% is too re- 
strictive and can miss some true BL Lac objects (Marcha 
et al. 1996, Laurent-Muehleisen et al. 1998, Rector et al. 
2000). In principle, this criterion may introduce in the 
sample some "optically passive" galaxies for which A is 
low due, for instance, to a particular evolutionary stage of 
the galaxy rather than to the presence of a BL Lac nu- 
cleus. Nevertheless, we have estimated (see Paper I) that 
the percentage of "normal" galaxies in the sample of BL 
Lacs defined according to the above definitions should be 
less than 10%, which means less than 5 spurious objects 
in the sample. A VLA survey of the BL Lacs and "opti- 
cally passive" galaxies discovered in the REX survey is in 
progress and should provide a compelling tool to separate 
BL Lac nuclei from "normal" galaxies. The redshift dis- 
tribution of the 37 BL Lac objects for which a redshift has 
been measured is presented in Figure 1. 

By assuming that the percentage of BL Lacs among 
the unidentified objects is the same as for the identified 
sources (i.e. 25%) we expect that only 5 BL Lacs will 
be found after the completion of the identification pro- 
cess. This number is in complete agreement with the one 
computed by comparing the overall spectral energy dis- 
tribution of the unidentified objects with that of the 55 
BL Lacs. Through this comparison we predict only 3-4 
missing BL Lacs. Given the small number of the expected 
missing objects, the conclusions of the analysis presented 
in this paper are not likely to change after the completion 
of the identification process. 
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Fig. 1. — The redshift distribution of the 37 BL Lac objects in 
the XB-REX sample with a spectroscopic measured redshift. 



1 http ; / / www. ast . cam .ac.uk/~apmcat/ 

^We have assumed O - E = 1.4 which is the mean value found in the sample 
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3 The Cosmological evolution: Results 

The high identification level of the XB-REX sample al- 
lows us to apply the K/Ki test extended to the multi- 
frequency selected samples (Avni & Bachall 1980). This 
test checks the imiformity of the spatial distribution of the 
objects in a flux-limited sample. A departure from a uni- 
form distribution is interpreted as due to an evolutionary 
behavior of the population under study. For each object, 
the two volumes Ve and Va are derived on the basis of the 
observed redshift (z) and the maximum redshift (zmax) at 
which the object would produce a flux equal to the flux 
limit of the survey. Since the XB-REX is characterized by 
3 flux limits, the most stringent value of z^ax is consid- 
ered, i.e.: 

ZmoK = {^maw> '^max'> '^max) 

where z^^^^, z^^^, '^max ^he maximum rcdshifts ob- 
tained respectively in the radio, optical and X-ray band. 
If the Ve/Va values are not uniformly distributed between 
and 1, there are indications for cither an evolutionary be- 
havior of the population under study or for the presence of 
selection effects biasing the sample. In particular the first 
requirement for a uniform distribution is that the average 
of the Ve/Va values is consistent with 0.5. The effective 
uniformity of the distribution must then be checked with, 
for instance, a Kolmogorov-Smirnov test. By using the 
sample of 55 BL Lac, we have computed the Ve/Va values 
for each source. For the objects without a redshift, we have 
assumed z=0.27 (which is the mean redshift of the sample). 
We have also assumed an = and ax = 1 for the spectral 
indices in the radio and in the X-ray band, respectively. In 
the optical band, BL Lacs show a wide range of spectral 
indices, with radio selected BL Lacs being steeper than X- 
ray selected BL Lacs (Falomo, Scarpa & Bersanelli 1994). 
Since the majority of our sources are more similar to the 
X-ray selected objects, we have used the mean value found 
by Falomo, Scarpa & Bersanelli (1994) for these sources 
{ao = 0.65). 

The resulting value of the < Ve/Va >, reported in Tab. 1 
(< Ve/Va >=0.484±0.039), is fully consistent with 0.5. 
The K-S test shows that the Ve/Vi values are consistent 
with a uniform distribution between and 1 (probabil- 
ity=76.3%). 

In order to check the stability of the results, we have 
performed the analysis under different assumptions. If we 
assume, for all the objects without a redshift determina- 
tion, the values taken from the extremes of the observed 
redshift distribution (Fig 1), i.e. 0.1 and 0.7, we obtain 
< Ve/Va > = 0.47 and 0.50 (±0.039) respectively. Thus, 
the result is not significantly affected by the number of 
objects without a redshift determination. 

The Ve/Va analysis may also be affected by a wrong as- 
sumption on the spectral slope. The expression used for 
the K-correction {K = {l + z)^"~^^), based on the assump- 
tion of a power-law spectrum with spectral index a, can 
mimic a luminosity evolution (or cancel a real one) If the 
value of a is over or underestimated. By assuming extreme 



values for the a in the 3 bands, namely Q!x=0.5, afi=- 
0.5, q;o=0.0 and ajc=2, q;_r=1, ao=2, the < Vg/Va > 
values range from 0.47 to 0.52. We conclude that the K- 
correction does not play a fundamental role in the analysis 
unless very extreme (and unrealistic) values of a are used. 

TABLE 1 

Results of the Vg/Va analysis 



Sample (N. objects) 


< Ve/Va > 


Error 


s2 


KB prob.(%)3 


Total (55) 


0.484 


0.039 


0.4 


76.3 


"classical" BL (46) 


0.461 


0.043 


0.9 


42.6 


aox <0.91 (22) 


0.488 


0.061 


0.2 


82.2 


anx <0.62 (25) 


0.469 


0.058 


0.5 


56.2 


EL AGNs (95) 


0.630 


0.030 


4.3 


«0.1 


simulated sample (64) 


0.423 


0.036 


2.1 


9.7 



1 Error = l/^/luN); 

The significativity defined as the distance from 0.5 in units of cr; 
^ The KS probability (null hypothesis = the Ve/Va distribution is 
uniform bewteen and 1). 



We have then investigated whether different types of 
objects show different cosmological behaviors. 

"Classical" BL Lacs. We have computed the < 
Ve/Va > by using only the "classical" BL Lacs, i.e. defined 
with the most restrictive limit on A (25%). The 46 "classi- 
cal" BL Lacs have a, < Ve/Va > lower than the entire sam- 
ple but still consistent with 0.5 (< Ve/Va >=0.461±0.043). 

HBLs We have considered only the 22 BL Lacs with 
an "extreme" X-ray/optical ratio {aox <0.91) typical of 
HBLs. We have also divided the sample according to the 
ratio between the X-ray and the radio flux and considered 
only those with anx < 0.62. In both cases the < Vg/Va > 
value is consistent with 0.5 (< Ve/Va >=0.488±0.061 and 
0.469±0.058 respectively). 

Emission Line AGNs. For comparison, we have com- 
puted the Ve/Va for the Emission Line AGNs found in the 
XB-REX sample. By using the sample of 95 EL AGN, we 
have computed the < Ve/Va > with the following assump- 
tions: an = 0.4 (the mean value found for the emission 
line AGNs in the REX survey, Caccianiga et al. 2000), ax 
= 1 and ao = 0.65. The result is indicative of a strong 
positive evolution (< Ve/Va >=0.630±0.030). We have 
then estimated the evolution by assuming a pure luminos- 
ity evolution of the form: L{z) = L{0){1 + z)'' under the 
hypothesis that the values of k in the three bands (radio, 
optical and X-ray) are the same. The best flt value is 
k =3.0 with a la interval of (2.5, 3.3). The value of k 
is intermediate between that found for the X-ray selected 
(fc=2.92, Delia Ceca et al. 1994) and the optically selected 
(e.g. A:=3.49, Padovani 1993) radio-loud AGNs. Both val- 
ues are consistent within la with the one found in the 
XB-REX sample. 
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4 The simulations 

The absence of cosmological evolution in BL Lacs is po- 
tentially at odds with the results presented in the litera- 
ture, claiming the existence of negative evolution (Morris 
et al. 1991; Woltcr ct al. 1994; Bade ct al. 1998; Giommi, 
Menna & Padovani 1999; Rector et al. 2000) or a posi- 
tive evolution (Stickel et al. 1991; Rector & Stocke 2001). 
Assuming a luminosity evolution of the form: 

L{z) = L{0)exp{cT) (1) 

where r is the look-back time, our result is consistent, at 
la confidence level, with a value of c between +2.0 and - 
3.5. Thus, both the negative (c=-7.0, Wolter et al. 1994) 
and positive (c=3.1, Stickel et al. 1991) best-fit parame- 
ters found in the literature fall outside the la range. 

It can be argued that the negative evolution of BL Lacs 
is not visible in the XB-REX sample simply because this 
sample is not deep enough. The flux limit of 4x10"^'^ erg 
s~^ cm~^ in the 0.5-2.0 keV band corresponds to a flux 
limit of 7.2x10-13 grg s'^ cm'^ in the 0.3-3.5 keV band 
(which is the EMSS energy band) assuming ax=l- This 
flux is higher than the lowest limit reached by the EMSS 
complete sample (5xl0~"'^'' erg cm^^) used by Mor- 
ris et al. (1991)3. The flux limit of the XB-REX is thus 
in the region where the BL Lac LogN-LogS measured in 
EMSS survey (Wolter et al. 1991) starts to flatten. Thus, 
the negative evolution found in the EMSS sample might 
simply not be observable in the XB-REX sample. 

In order to verify this hypothesis we have analyzed a 
simulated sample (see Paper I) created under the assump- 
tion of negative evolution of the exponential form given 
by Eq. 1 with the best-fit value of c = —7.0 found by 
Wolter et al. (1994). If we apply the Vg/Va analysis to 
the entire simulated sample and we compute the evolu- 
tionary parameter, we find a best-fit value of —6.8, i.e. in 
excellent agreement with the input value. This first test 
demonstrates the validity of the method when applied to 
a multifrequency survey. 

If we apply to the simulated sample the same radio, 
optical and X-ray flux limits used for the XB-REX, the 
resulting sample contains 65 objects. This number is in 
excellent agreement with the observed number of objects 
if we consider that 10% of the objects are expected to 
bo missed in the positional radio/X-ray correlation as ox- 
plained in Paper I and an additional 9% is expected to 
be missed due to the identiflcation level (91%) achieved so 
far in the XB-REX sample. Considering those corrections, 
the number of BL Lacs expected in the XB-REX sample 
is about 0.9x0.91x65 = 53 that should be compared with 
the observed number of 55. 

The < Ve/Va > analysis of the simulated sample shows 
the presence of the negative evolution (< Ve/Va >= 0.423 
± 0.036, SCO Tab. 1). The Ve/Va distribution is not luiiform 
although the KS test rejects the hypothesis of no-evolution 
only marginally (probability of 9.7%). By assuming the 

^the flux limit of 5x10"^^ erg cm~^ is actually reached only 
limit corresponding to the 90% of the area is much higher (~2.5xl0~ 



same form of the evolution used as input for the simula- 
tions, we obtain a best-fit value (c=-4.5 with la range 
of -7.2, -2.2) which is lower than the input value (-7.0) 
although consistent at Icr. Thus, the negative evolution 
is recovered but it is weaker than the one given as input. 
This result is probably due to the fact that, as the flux 
limit increases, the negative evolution becomes less and 
less evident since the euclidean part of the LogN-LogS is 
sampled. 

In conclusion our simulations show that, at the X-ray 
flux limit of the XB-REX sample, the negative evolution 
seen in the EMSS sample should be less extreme but still 
detectable. The result of the analysis presented here may 
indicate that, if present, the negative evolution is weaker 
than that estimated from the EMSS sample. 

5 Discussion and conclusions 

We have presented the result of the Ve/Va test on a 
sample of 55 BL Lacs extracted from the REX survey. No 
evidence of positive or negative evolution is found. The 
analysis presented in this paper is rather stable, i.e. it is 
not affected much by the fact that some sources do not 
have a rcdshift determination (indeed, the Ve/Va test does 
not depend very much on the redshift of the sources). 

Simulations based on the EMSS sample has demon- 
strated that, at the X-ray flux limit of the XB-REX sam- 
ple, the negative evolution found by Wolter et al. (1994) 
is expected to be less evident although still detectable. 
More puzzling is the discrepancy between the results found 
in the XB-REX sample and the strong negative evolution 
presented in Bade et al. (1998). Since the RASS sample 
has a higher X-ray flux limit than the XB-REX sample, 
the strong negative evolution found in the RASS sample 
should have been clearly detected in our analysis. 

A possible explanation of this apparent contradiction 
is that the evolutionary behavior is not the same for all 
the different sub-classes of BL Lacs. As suggested by dif- 
ferent authors, there may be a continuous trend of the 
evolutionary properties within the BL Lac class, with the 
most extreme HBLs evolving negatively, from one side, 
and the most extreme LBLs evolving positively, from the 
other side. Depending on the flux limit (s) of a given sur- 
vey and/or to its selection criteria, a sample may contain 
a different percentage of HBLs and LBLs and this can pro- 
duce different results in the Ve/Va analysis. In Figure 2 
the anx distribution of the 55 BL Lacs selected in the XB- 
REX sample (tick line) is presented and compared to the 
same distribution for the EMSS (shaded histogram on the 
left side), the 1 Jy (filled histogram on the right side) and 
the RASS (dashed line) BL Lac samples. It is clear that 
the XB-REX BL Lacs are covering a range of anx which 
is similar to that observed in the X-ray selected samples 
(EMSS and RASS) although it does not reach the very 
low end of the distribution {aux <0.5). Hence, we may 

a fraction (58%) of the area covered by the EMSS survey. The flux 
erg s~l cm~^) 
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expect qualitatively that the evolutionary behavior of the 
XB-REX BL Lacs is similar to that observed in the X- 
ray selected samples but less extreme since the very low 
o:rx objects are not included in the sample. Since we do 
not observe evidence for evolution even when we consider 
the objects with a^x <0.62, the conclusion is that the 
extreme negative evolution must be carried out mainly by 
the minority of the objects with urx most extreme, not se- 
lected in our sample. Surveys more specifically designed to 
find the extreme HBLs, like the Sedentary multifrequcncy 
survey, should be in the best position to observe the signs 
of a negative evolution (Giommi et al. 2001). 
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Fig. 2. — The radio-to-X-ray spectral index (ajjx) distribution 

of the BL Lac objects in different samples: XB-REX (tick line), 
EMSS (shaded histogram on the left side), IJy (filled histogram 
on the right side) and RASS (dashed line) 

In any case, the analysis presented in this paper con- 
firms that BL Lac objects do not follow the same evo- 
lutionary behavior of the emission line AGNs. Possible 
explanations for this difference have been recently pro- 
posed by Cavaliere & Malquori (1999) on the basis of a 
different accretion rate on the central massive black-hole: 



BL Lac objects may be characterized by a low accretion 
rate {m = Mc^/Le ~10~^) that produces a very mild 
luminosity evolution in contrast with what happens in the 
Flat Spectrum Radio Quasars (FSRQ), where the accre- 
tion is close to the Eddington limit (m ~1) producing 
a much faster luminosity evolution. The predicted time 
scale for the luminosity evolution in the BL Lac popu- 
lation is expected to give a < Ve/Va > consistent with 
0.5, as observed in the XB-REX sample. Superimposed 
to this effect, a possible transition from a FSRQ into a 
BL Lac object, due to a reduction of the accretion rate 
at low redshifts, may cause negative (density) evolution 
in BL Lac objects (objects more numerous now than at 
high redshift) as described in Cavaliere & D'Elia (2001). 
A similar evolutionary link between FSRQ and BL Lac has 
been proposed by Bottchcr & Dcrmcr (2001) where the re- 
duction of the black-hole accretion power with time may 
switch a FSRQ into a BL Lac object following the sequence 
FSRQ^LBL^HBL. As discussed before, the analysis of 
the XB-REX sample indicates that only a mild negative 
evolution is consistent with the data although only a nar- 
row range of BL Lac "flavors" is actually covered in this 
sample (see Fig. 2). After the completion of identification 
process of the entire REX survey, which contains 7 times 
more sources than the XB-REX sub-sample, a wider range 
of BL Lac types will be covered allowing a more complete 
study of the evolutionary behaviors within the BL Lac 
class. 
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TABLE 1 
The XB-REX sample of BL Lacs 
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" This object has been presented in Wolter et al. (1997) with the name REXJ1231. 4+1421. 
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TABLE 1 
(continued) 
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Column 1: Name, based on the radio position; 

Column 2: The unabsorbed X-ray flux in the 0.5-2.0 keV energy band; 

Column 3: The integrated radio flux density at 1.4 GHz (from the NVSS catalog); 

Column 4: The optical blue magnitude. We have used preferentially the APM O magnitude. 

When not available, a blue magnitude has been estimated in the following ways: 

(^) using the APM E magnitude, as described in the text 

(^) using the V magnitude published in Padovani & Giommi (1995) assuming B-V = 0.6 (mean value for the BL Lacs). 
(^) using the B magnitude published in Nass et al. (1996) 

Column 5: The redshift. Tentative redshifts are indicated with (*). The redshifts are either computed on the basis of our spectra or they are taken from 
NED. 

(1) Bauer et al. (2000) reports z=0.2106 

(2) Cao, Wei & Hu (1999) reports z=0.27 

(3) redshift taken from Padovani & Giommi (1995). Difi'erent redshift in NED (z=0.13) 

(*) redshift taken from Padovani & Giommi (1995). Different redshift in NED (z=0.182) taken from Bade et al. (1996); 
(5) Bauer et al. (2000) reports z=0.469 

Column 6: A flag indicating that the object has been observed spectroscopically by us (flag="Y"). 

The flag "Y" does not mean necessarily that the object is newly discovered. If a classification as BL Lac is already present in the literature an alternative 
name is reported in column 7. 

Column 7: An alternative name from the literature, when available. 

Column 8: type of BL Lac: "BLC"= BL Lac candidate, i.e. the value of A is between 25% and 40%. 
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